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' 1. Introduction 
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We present combined results on the search for neutral supersymmetric Higgs bosons with data collected at the 
CDF and DO experiments. Data were collected in proton-antiproton collisions at y's = 1.96 TeV during Run II 
of the Tevatron. The searches considered cover the main production and decay mechanisms for Higgs bosons in 
tau lepton and bottom quark final states. The resulting combination is interpreted in the context of different 
scenarios within the Minimal Supersymmetric Standard Model. 



Supersymmetry (SUSY) is an extension to the standard model (SM) that predicts an additional symmetry 
^ ^ between bosons and fermions. This results in several advantages over the SM, such as the introduction of a 
^ i dark matter candidate, a solution to the hierarchy problem and a potential for GUT scale unification. The 
i' Minimal Supersymmetric Standard Model, MSSM, [l| is the simplest version of a SUSY theory. It predicts, 
1^ . after symmetry breaking, five Higgs bosons of which three are neutral, these are denoted h, H, A, H~ . At 
i ' tree level the Higgs sector of the MSSM can be described by two parameters, which are chosen to be tan/3, 
the ratio of the vacuum expectation values of the two Higgs doublets, and Ma, the mass of the pseudo-scalar 
—i Higgs boson. In the MSSM the coupling of the neutral Higgs to down-type quarks and leptons is enhanced by a 
^ factor of tan/3 and the corresponding coupling to the up-type quarks and leptons is suppressed. This means that 
■ decay modes to bottom quarks and tau leptons are of specific interest, with a predicted decay rate of around 
' 90% to bottom quarks and 10% to tau leptons. At large values of tan/3 two of the three neutral Higgs bosons 
I , are approximately degenerate in mass which leads to an approximate doubling of the production cross section. 
• ' Beyond tree level, radiative corrections bring in dependance on more than just Ma and tan/3. Therefore, limits 
] are set for specific benchmark regimes in the Ma, tan/3 plane. For regions of the parameter space with low 
values of Ma and large values of tan/3 the Tevatron can set strong limits that complement the previous searches 
by the LEP experiments 01 ■ 



2. Combined CDF and DO MSSM combination 

The combined results in the di-tau final state from both DO and CDF are presented using up to 2.2 fb^^ and 
1.8 fb^^ of integrated luminosity, respectively. The searches are described in ^SnSi] and detailed information on 
the DO and CDF detectors is available in [1, 0]- The CDF analysis studied three final states, where either one 
tau decays leptonically to a electron or muon and one hadronically or both decay leptonically, TeThad,T^Thad 
and Tf.T^. At DO a 1.0 fb^^ analysis of these three final states was combined with an additional 1.2 fb""'^ analysis 
that studied just the TnThad final state. More details on this combination are available in Q. 



2.1. Analysis Summary 

Electrons are identified through their characteristic energy deposits in the calorimeter and are required to 
be isolated and have a reconstructed track matching this energy deposit deposition in the calorimeter. For 
reconstructed muon objects, hits in the muon chambers are matched to tracks in the central tracking chambers. 
These are then required to be isolated in both the calorimeter and the central tracking detectors. DO and CDF 
have different tools to identify hadronic decays of taus. At CDF narrow isolated clusters in the calorimeter 
with either one or three tracks are identified and a variable cone algorithm is used to reconstruct these with 
strict isolation requirements in order to suppress quark and gluon jets. Hadronic tau decays at DO are first 
split into three types, based on the energy clusters in the calorimeters and the number of tracks. Type 1 and 2 
are 1-prong decays with energy deposited in the hadronic calorimeter (type 1) or in both the electromagnetic 
and hadronic calorimeters (type 2). Type 3 taus are 3-prong decays with an invariant mass below 1.7 GeV and 
energy deposits in the calorimeters. For each of the different types a Neural Network, NN, is trained, using 
Z — > TT decays as signal and multijet events predicted from data as background, to separate hadronic tau 
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Table I: Expected numbers ol background and signal events for rr final states. 



CDF L = 1.8 fb^^ 




DO L < 2.2 fb-i 


'?~e '^h a d 'Tm u '^h a d 


Z ^TT 

Z ee/Z — ^ fifi 
diboson+tt 
multijet+W — Iv 


605 ± 51 1378 ± 117 1353 ± 116 

19.4 ± 5.7 70 ± 10 107 ± 13 

20.5 ± 7.0 8.2 ± 4.2 6.6 ± 3.7 
57.1 ± 13.5 467 ± 73 285 ± 46 




212 ± 20 581 ± 5 2153 ± 156 
10.4 ±1.3 31 ± 2 66 ± 8 
6.1 ± 0.6 3.1 ± 0.3 16 ± 3 
37.9 ± 7.7 374 ± 48 216 ± 41 


Total Background 


702 ± 55 1922 ± 141 1752 ± 129 




266 ± 22 989 ± 82 2451 ± 162 


Data 


726 1979 1666 




274 1034 2340 



decays from jets. An additional NN trained on electron Monte Carlo events is used to remove backgrounds from 
electrons faking type 2 taus. 

Monte Carlo simulations are used to determine the acceptance of the signal. The PYTHIA Q event generator 
is used with CTEQ5L parton sets [H at CDF and CTEQ6L at DO. TAUOLA [ll| is used to simulate the 
final state tau lepton decays. For CDF both gg and hb production modes are considered, DO just considers 
gg production. The acceptances are very similar for both production modes. The response of the detector 
is modeled by GEANT [12|, the SM backgrounds are generated by PYTHIA, except tt which is generated 
by COMPHEP with PYTHIA and W/Z samples with additional jets in the sample which is generated 
by ALPGEN with PYTHIA used as normalization. Diboson and tt are normalized to next-to-leading 
order and Z/^* samples are normalized to next-to-next-to- leading order [l3]- DO uses a inclusive electron 

and muon trigger and CDF a lepton plus track trigger. After reconstruction, candidate events must have two 
isolated leptons of opposite charge. The main backgrounds are Z tt, fj^fj,, ee, multijet, W ev, ^u,tv, 
diboson and it. For the TeThad and T^Thad channels the following cuts are applied at CDF; the muon or electron 
is required to be isolated with a pr cut of > 10 GeV, the tau is required to have px > 15 GeV for a one-prong 
tau and > 20 GeV for three-prong taus. At DO the muon or electron is required to be isolated with a pt cut 
of > 15 GeV , the tau is required to have pr > 16.5 GeV for a one-prong tau and > 22 GeV for a three-prong 
tau. Additional cuts are placed on the relative directions of the tau, the I^t and the transverse mass. At CDF 
a cut is also placed on the scalar sum of the transverse momenta, defined as Ht ~ \Pt^\ + \Pt""'\ + Wt\ > 55 
GeV. For the TsT^ analysis, events with a central muon and electron are selected, with mm{E^,pt^) > 6 GeV, 
max{E^,pt^) > 10 GeV and \E^\ + \p!^\ > 30 GeV. At DO it is required that > lOGeV, p^ > 12 GeV, 
invariant mass of the e/z pair > 20 GeV and -I- -I- \1^t\ > 65 GeV. The expected yields for signal and 
background events can be seen in Table I along with the observed data events. 

2.2. Combination 

Two different methods are used to set the limits to be confident that the final results do not depend on 
the details of the statistical formulation. These are the Bayesian and Modified Frequentist approaches and 
rely not just on the integrated yields of the analyses but on the shapes of the final discriminants. Systematic 
uncertainties appear as both uncertainties on the number of expected signal and background and distributions 
of the final discriminates. For both DO and CDF the visible mass defined as the invariant mass of the visible 
tau decay products and I^t was used as the final discriminant. 

2.3. Systematic Uncertainties 

The systematic uncertainties that are correlated between the two experiments include the uncertainties on 
the integrated luminosity, on the rate for tt production, signal and diboson production. As both CDF and DO 
extract the multijet background from the data using two different techniques, the uncertainties on the multijet 
prediction is uncorrelated between the two experiments. Similarly the methods of calibrating the fake lepton 
rates, b-tag efficiencies and mistag rates are performed separately by each experiment and are uncorrelated. 

2.4. Combined Results 

The two methods used to set the limits show good agreement with each other with deviations at less than 
the 10% level. These results can be seen in Fig. 1 using the CLg method over the range 90 < Ma < 200 GeV. 
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The observed limits are in good agreement with the expectation with no significant excess observed. 




Tevatron Runll Preliminary 1.0-2.2fb'^ 



,n-1 

^"lOO 110 120 130 140 150 160 170 180 190 200 

[GeV] 

Figure 1: 95% confidence limits on the cross section x branching ratio of a neutral MSSM Higgs boson. The solid and 
dashed lines show the observed and expected limits respectively. The yellow and blue bands show the la and 2a bands. 

These limits can be interpreted within the MSSM in specific benchmark scenarios. Four specific scenarios are 
defined by setting the values of Msusy, the mass scale of the quarks, /i, the Higgs sector bilinear coupling, M2, 
the gaugino mass term, At, the trilinear coupling of the stop sector, Ab, the trilinear coupling of the sbottom 
sector and rUg, the gluino mass term. The maximal mixing scenario m™"^ is defined by setting 

Msusy = 1 TeV, a* = 200 GeV, M2 = 200 GeV 

Xt = 2MsusY, Ab=At, rug^ O.SMsusy- (1) 

The no-mixing scenario is defined by requiring Msusy — 2 TeV and Xt = 0, with the rest of the parameters 
unchanged from the m™"^ scenario. The two additional scenarios are created be changing the sign of the n 
parameter. The 95% confidence limits set in these scenarios can be seen in Figure 2 for the maximal and no 
mixing scenarios and a positive /i parameter. FEYNHIGGS [l^ was used to calculate the signal cross section 
and the branching fraction with no theoretical uncertainties considered. For large values of tan/3 the width of 
the Higgs boson can be compared to the experimental resolution. This width dependance is not taken into 
account when calculating the limits but in the region where the limits are set it is not expected to be a large 
effect [H. 



3. DO MSSM combination 

The combined results from DO on the search of a neutral Higgs boson produced in association with a h quark, 
in the bh — )> hrr and bh bbb channels, using an integrated luminosity of 7.3 fb~^ and 5.2 fb^-'^, respectively, is 
presented. Limits are set at the 95% confidence level on the production cross section of a neutral Higgs boson 
produced with one or more b quarks, (j{gb hb{b)), for Higgs boson masses 90 < Ma < 300 GeV, where it is 
assumed that BR{h bb) + BR{h tt) — 1 and there is a narrow Higgs width. This combination is described 
in detail in [l^] • 

3.1. Analysis Summary 

Detailed descriptions of the two combined analyses bh — > brr and bh bbb are available in [2l|, [l^. The 
muon and hadronic tan identification is performed as described in Section 2.1. For an object to be selected as a 
6-jet candidate it must pass as set of quality criteria. Within a cone of Ai? — ^ (Ar/)^ + (Ac/))^ < 0.5 around the 
jet axis it must have two or more tracks. A Neural Network algorithm trained on lifetime information from the 
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Figure 2: 95% confidence limits in tiie Ma tan/3 plane, for tlie two benclimark scenarios and /i > 0. Tlie black line 
denotes the observed limit, the grey line the expected limit and the hatched yellow and blue regions the Icr and 2a bands. 
The green region is the region excluded by LEP. 



track input parameters and secondary vertices is then used to further select events [231] . For more information 
on jet reconstruction and energy scale determination, see [2^ . 

To determine the signal acceptance, the gh hb{b) process is simulated using PYTHIA in the 5-flavor 
scheme at leading order. NLO corrections are calculated using MCFM 25] and are applied as weights on the 
transverse momentum, pr, and the r/ distributions. TAUOLA is used to simulate the final state tau lepton 
decays and GEANT to model the detector response. Weights are applied to correct for the differences between 
the actual detector response and the simulation. All backgrounds are simulated by ALPGEN, with PYTHIA 
used to model the hadronization and showering, apart from the diboson processes which are simulated by 
PYTHIA only. The diboson and tt samples are normalized to next-to-leading order, whereas Z/'y* samples are 
normalized to next-to-next-to-leading order. 



3.1.1. bh bTT 

For the bh brr analysis, where one tau is required to decay hadronically and one to a muon, an inclusive 
trigger approach is used which is a combination of the signal muon, electron, jet, tau, muon -I- jet and muon -|- 
tau triggers. One isolated muon is required with a matching central track, which must pass the cuts: px > 15 
GeV, \ridet \ < 1.6 . Events with greater than one muon are rejected to suppress the Z — > fifj, background. The 
hadronic tau candidates must be isolated with the cuts: pt > 15 GeV, |?77-| < 2.5. A cut is also applied on the 
tau Neural Network output, NNr, which selects hadronic tau candidates with a efficiency of 65% while rejecting 
jets with a efficiency of around 99%. Events must also have a 6-jet which passes the following criteria; isolated 
from the tau and muon, px > 15 GeV, \ri\ < 2.5, \r]det \ < 2.5. In addition a cut is made on the b tagging Neural 
Network, NN;,, which selects 6-jet candidates with an efficiency of around 65% and a fake rate of 5%. The 
main backgrounds are multijet, ti and Z ^ tt produced with heavy flavor jets, with small contributions from 
VF-|-jets, Z TT plus light jets, Z -> II, single top production and diboson production. For VF-|-jets the Monte 
Carlo is normalized to data using a control sample. Multijet background is estimated from data using inverted 
cuts to select a multijet-rich sample. To further reject the main backgrounds, multijet and tt. A discriminant 
is trained on the multijet background, DmJi and a NN on the tt background, Dj^. The final discriminant that 
is used for the limit setting is a likelihood discriminant, Df, which is constructed from Dmj and Dii as well 
as the ^-tagging NN, and a variable relating to the decay kinematics, Mhat- The distribution for Df\s shown 
in Figure 3 for all tau types and a 110 GeV Higgs signal (for the statistical analysis the tau types are treated 
separately) . 

3.1.2. hh bbb 

For the bh — > bbb analysis a trigger specifically designed to select events with at least 3 jets is used. This 
trigger is determined to be around 60% efficient for a 150 GeV Higgs signal when measured with respect to 
events with 3 or 4 reconstructed jets. Events with 3 or 4 jets are selected, 3 of which must satisfy pr > 20 GeV, 
jyl < 2.4 criteria and must pass tight NN b tagging cuts. The efficiency of the b tagging NN is around 50%, 
with a fake rate for light jets of around 0.5-1.5%. An additional cut on the pr of the two leading jets of greater 
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D0, L = 7.3 fb-^ ^ data 




D, [m.^ = 110GeV/c2] 

Figure 3: The likeihood discriminant used in the bh — >■ 6rr analysis, for a 110 GeV Higgs mass. The stacked histograms 
show the contribution from the different backgrounds and the hatched histogram show the predicted Higgs signal. 

than 20 GeV is also applied. The background is dominated by heavy flavor jets with a smaller contribution 
coming from light flavor jets. This multijet background is predicted using a data driven method. To enhance 
the selection of the signal over the background, a likelihood discriminant is used based on a set of kinematic 
variables for which the background is well modeled. Two likelihood discriminants are trained, one for the low 
mass region 90 < Ma < 140 GeV and one for the high Higgs mass region 140 < Ma < 260 GeV. The jet pair 
with the highest likelihood is selected and used to construct an invariant mass which is used as the input to the 
statistical analysis. This distribution can be seen in Figure 4. 

3.2. Systematic Uncertainties 

Systematic uncertainties are either included as constant systematics on the yield of signal or background 
events or as a shape systematics, the value of which changes across the different bins of the distribution. The 
dominant sources of systematic uncertainties for bh brr analysis are as follows. For the Z+jets samples the 
uncertainness are estimated from Z — > fifi events, these include the Z+jets normilazation uncertainties, the 
inclusive trigger efficiency (applied to all background samples) and a shape uncertainty from the modeling of 
the Z boson kinematics. The other background samples are dominated by the uncertainties on, the hadronic 
tau reconstruction efficiency, the luminosity determination, the muon reconstruction efficiency, the single muon 
trigger efficiency, the ti and diboson production cross sections. The determination on the multijet background 
has a large uncertainty of 10-40%. Additional dominant shape uncertainties come from the jet energy scale 
determination and the 6-tagging efficiency modeling. Apart from the uncertainties on the multijet production 
and the hadronic tau reconstruction efficiency all these uncertainties are taken to be 100% correlated across the 
tau types. 

In the case of bh — )■ bbb only shape systematics are considered. The dominant background arises from the 
measurement of the rate at which light partons fake heavy flavor jets and the b tagging efficiency. In contrast, for 
the signal the dominant sources of uncertainness are the jet energy scale and the b tagging efficiency. Between 
the two analyses only the b tagging efficiency and the jet modeling systematics are correlated and are taken to 
be 100% correlated. 

3.3. Combined Results 

Limits are set using the modified frequentist or CL5 technique [26| , with more details on this available in (20| . 
This is done both within the the framework of the MSSM at the benchmark points as described in Section 2.4 
and in a less model independent way on the cross section x branching ratio. This is not normally possible when 
combining channels with different production modes or final states, as the relative signal yields in each channel 
will generally depend on the specific model being studied. But as the two channels share a production mode, 
and the Higgs width is narrower than the experiment resolution for the majority of the interesting phase space, 
model independent limits can be derived. The only additional assumption that must be made is that the Higgs 



6 



Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011 




D0 , 5.2 fb-' 
a) 3 jet 

Low-mass likelihood 
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Figure 4: The invariant mass of the highest likelihood pair of b quarks. The line shows the predicted background, with 
the heavy flavor contribution in blue. The crosses show the data with the difference between the data and the expected 
background shown below. 



boson only decays to tau leptons and h quarks. This model independent limit can be seen in Figure 5, for 3 
choices of the tau branching ratio, BR = 6%, 10%, 14%. For masses up to around 180 GeV, the bh — > hrr 
channel tends to dominate, beyond this value the bh — > bbb becomes increasingly important as can be seem from 
the decrease in the dependence of the limits on the tau branching fraction at higher values. 




Figure 5: Combined 95% confidence rates on the cross section. The observed (solid) and expected (dashed) are shown 
for three assumptions of the tau branching ratio. The shaded bands show the theoretical uncertainty an error of 15% 
for the cross section for the m^""^ scenario for two different assumptions of and tan/J = 40. 

In Figure 6 the limits are set in the mj^"^ and no mixing scenarios as described in Section 2.4, for a positive 
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value of tan/?. The experimental uncertainties an additional 15% uncertainly on the production cross section is 
included. This has contributions from the scale variations and the PDF uncertainties. The SM cross section of 
the gh — >■ hh process as taken from MCFM where FEYNHIGGS is used to calculate the branching fraction so 
that the SM cross section can be corrected to the appropriate cr x BR. For large values of tan/? the width of 
the Higgs boson is comparable to the experiment resolution, this effect is taken into account using the method 
described in j^^]. These limits are the best limits so far produced at the Tevatron. 




Figure 6: 95% confidence limits in the M^-tan/? plane, for the two benchmark scenarios and a positive jj, value. 
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